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Exact ground states, charge densities and excitation energies are found using valence bond meth- 
ods for A'^-site modified Hubbard models with uniform spacing. At the neutral-ionic transition 
(NIT), the ground state has a symmetry crossover in 4n, 4n-|-2 rings with periodic and antiperiodic 
boundary conditions, respectively. Large site energies A stabilize a paired state of the half-filled 
f^ , chain, while large U stabilizes a covalent state. Finite transfer integrals t shift the NIT to the cova- 

lent side oi U — 2A. Exact results to A*' = 16 in the full basis and to A'^ = 22 in a restricted basis 
for large U, A are extrapolated to obtain the crossover and charge density of extended chains. The 
flj . modified Hubbard model has a continuous NIT between a diamagnetic band insulator on the paired 

side and a paramagnetic Mott insulator on the covalent side. The singlet-triplet (ST) , singlet-singlet 
(SS) and charge gaps for finite N indicate that the ST and SS gaps close at the NIT with increasing 
U and that the charge gap vanishes only there. Finite-A*' excitations constrain all singularities to 
± O.lt of the symmetry crossover. The NIT is interpreted as a localized ground state (GS) with 
finite gaps on the paired side and an extended GS with vanishing ST and SS gaps on the covalent 
side. The charge gap and charge stiffness indicate a metallic GS at the transition that, however, 
1-^ I is unconditionally unstable to dimerization. Finite A breaks electron-hole (e-h) symmetry, but the 

f— ! ■ modified Hubbard model has an extended e-h symmetry and strong mixing of spin and charge exci- 

Q ' tations is limited to a few t's about the NIT. Exact finite-size results complement other approaches 

O , to valence or ferroelectric transitions in organic charge-transfer salts or in inorganic oxides, and to 

electron-vibration coupling and structural instabilities in one-dimensional systems. 
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I. INTRODUCTION 

McConncU and coworkersO explained the sharp sep- 
aration of organic charge transfer (CT) complexes into 
diamagnetic and paramagnetic by proposing that weak 
TT-donors (D) and acceptors (A) form neutral complexes 
of molecules, while strong donors and acceptors crystal- 
lize as ion radicals D+ and A~ . These planar conjugated 
systems form one-dimensional structures, either as mixed 
...DADA... stacks in CT complexes or as segregated 
stacks in ion-radical saltaj. The simplest approximation 
for a crossover between DA and D+A^ ground states is 
M = Ej—Ea, where M is the Madelung energy, Ej is the 
ionization potential of the donor and Ea is the electron 
affinity of the acceptor. Although M is inherently long 
ranged, the systems are quasi one-dimensional bjj-, virtue 
of TT-overlap restricted to stacks. Strebel and SoosB intro- 
duced the modified Hubbard model with transfer integral 
t = -{DAIHID+A-) for CT complexes and studied the 
crossover in the random phase approximation. Finite t 
leads to mixing and to partial ionicity D''+A''~ in the 
ground state (GS), with Ap < 1 at the crossover. The 
modified Hubbard model, Eq. (n|) below, has proved to 
be extremely rich and widely applicable. It describes any 
valence transition, is a special case of important solid- 
state models, and provides the starting point for electron- 
phonon (e-ph) coupling. Its scope is still growing and at- 
tracting new theoretical and computational approaches 
to the crossover region, M ^ Ej — Ea- We present in 
this paper exact solutions of finite-size systems, including 
low- lying excitations. 

The neutral-ionic transition (NIT) originates with t 
TTF-CA complex studied by Torrance and coworker; 
with D = tetrathiafulvalene and A = chloranil. TTF- 
CA is neutral at room temperature, with p ~ 0.3, and 
has a transition at T ~ 81K to an ionic state with p ^ 
0.7 that, moreover, is dimerized. The uniform TTF- 
CA spacing above 81K becomes alternating (...ii,i2---) 
in the ionic phase, and partial ionicity is determined 
spectroscopicallyQ. The structural change shows the fun- 
damental role of lattice phonons and the Peierls insta- 
bility of the paramagnetic phase. The alternating phase 
is potentially ferroelectric and the system may be metal- 
lic at the NIT. Such features are common to ferromag- 
netic oxides, and in this context they have recently been 
discussedS in terms of the modified Hubbard model. 
The interplay of electron-electron (c-e) and e-ph inter- 
actions can generate either continuous or discontinuous 
ionicity changes. Long-range Coulomb interactions can 
generateLTlJ discoatinuous p variations at the NIT as well 
as stropcly affecuJ the dimerization instability. 

Riceli3 pointed out the strong infrared activity of to- 
tally symmetric molecular vibrations through coupling 
to charge fluctuations. These on-site (Holstein) phonons 
also participate in the NIT. They condense at the tran- 
sition and produce discontinuous p variations above a 
critical coupling strengthtl. Electron-molecular-vibration 



coupling provides the basis for the spectroscopic deter- 
mination of the ionicity, p, of the D''+A''~ GS as well as 
of the local symmetry, making vibrational spectroscopy 
an useful ptool to follow charge and structural phase 
transitional^. Joint theoretical and expprimental anal- 
ysis of electronic and vibrational spectrallE-allowed for a 
systematic characterization of several saltstj. 

The modified Hubbard model adds site energies ±A to 
a Hubbard chain with uniform spacing, 

Ho{t, A,U) = -J2{t «.a.+i,. + a++l,,a^,.) - A(-l)'a+,a,,,} 

i,(T 

(1) 

D and A are at odd and even i, respectively, in the con- 
text of CT complexes, and i > 0, A > can be taken 
without loss of generality. The half-filled case, with one 
electron per site, is by far the most important. We con- 
sider Hq at this filling for uniform t and equal [/ > for 
donors and acceptors. The electron density on D sites is 
related to the GS energy. 
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The ground states of extended systems are not known 
exactly for A ^ 0. Approximate solutions beyond 
mean field have been .woposed along several liaes: ex- 
act diagonalizationD~lj'La, qiiantum Monte Carldlj, renpt- 
malization group methodalj, and continuum modelsliZI. 
There is broad agreement as well as open or disputed 
points mentioned below. 

Hubbard models are readily generalized and provide a 
unified approach to quantum cell models that need not be 
low-dimensional. In the context of Eq. (Il|) , we note that 
A can incorpo|uate the Madelung energy in a mean-field 
approximatioio or coupiJBS to Holstein phonons in the 
adiabatic approximatiodj'til. In either case, the effective 
A depends on the GS ionicity and the NIT becomes dis- 
continuous above a critical coupling. The model is then 
nonlinear and has wider applications to susceptibilitiesllS. 
At A = 0, if i is linearly expanded around the eqiii- 
librium bond-length, we get a Peierls-Hubbard modelEj, 
and a direct connection to models for ion-radical salts 
such as TTF-TCNQ with segregated stacks. Alternating 
transfer integrals i(l ± S) are found on the ionic sidmLin 
many segregated stacks and in conjugated polymersotiJ. 
Theoretical interest in Ho{t, A, U) and its variants lies in 
the interplay of e-e and e-ph interactions and their role 
in the sfa'pctural instabilities of 1-D materials, broadly 
definedllTB. 

The i = GS of the Hamiltonian in Eq. (|l|) is sketched 
in Fig. 1 to illustrate some basic features. The electrons 
are paired on D (odd) sites for A > U/2, paired on A 
(even) sites for A < —U/2, and singly occupy all sites 
in between; covalent states have n^ = 1 at all i and spin 
degeneracy of 2^. The valence transitions at U — ±2 A 
for C/ > transfer one electron in the GS and p = 2 — njj 



changes discontinuously. Since the paired GS in Fig. 1 
are nondegenerate, we expect finite gaps for spin, optical 
and charge carrying excitations. The covalent GS, on the 
other hand, has vanishing spin gaps. 

H(){t, A, U) with finite t is a one-dimensional metal at 
A = [/ = 0, a band insulator for A > 0, U = 0, and 
a Hubbard model for A = 0, ?7 > 0. Finite t leads to 
continuous no at the NIT when A is not a function of 
no. Although the ionicity is continuous, the NIT be- 
tween band and Mott insulators is a true quantum phase 
transition al|-.T=OK as signaled by the closing of tripletCp 
and singletaB gaps and by the unconditional instabilitylS't 
to dimerization on the covalent side. The nature of the 
transition between two insulators has revived interest in 
the modified Hubbard model in connection wf^ local- 
ization and conductivity in correlated systern^iL Strong 
charge fluctuations inducedj-hy lattice motionOo and re- 
lated structural instabilitiescj have been rediscovered and 
underlined. Finite t generates correlated states of the 
Hamiltonian in Eq. (||) that differ fundameutall^ ^^Pi U^ 
t = limit. In spite of sustained researchuijo ljEj lj, 
no definitive picture has emerged for the T=OK phase 
diagram of the simple model in Eq. {M . 

We present in this paper exact solutions ofpffo(t, A, U) 
for finite N using valence bond (VB) methodaSil that were 
originally developedQ for CT complexes. Total spin S is 
conserved in all versions of Eq. ffl). VB diagrams with 
specified pairing of sites with n^ = 1 form a large but 
complete basis for any N. The scope of finite-iV re- 
sults is decisively extended by using both periodic and 
antiperiodic boundary conditions and the symmetries of 
Hq. The oligomers in Section 2 reach N = 16 in the full 
basis or iV = 22 in a restricted basis without D^+ or A^^ 
sites, compared to TV ~ 10 in previous studies. Exact 
excitations near the NIT are related in Section 3 to the 
opening of gaps and interpreted as due to localization on 
the paired side. The GS is metallic at the NIT according 
to the charge gap and charge stiffness. The NIT marks 
the boundary between a localized GS for A >> U, and 
a delocalized GS with vanishing excitation energies as in 
Hubbard models at A = 0. We shall briefly mention the 
role of e-ph coupling and intersite e-e interactions, but 
defer detailed analysis to subsequent publications. 



II. SYMMETRY CROSSOVER AND CHARGE 
DENSITY 

We consider GS properties of Ho{t,A,U) at half fill- 
ing, one electron per site, and uniform < = 1. General 
solutions of the Hamiltonian in Eq. (|^) are restricted to 
finite N, where eigenstates and energies are accessible. 
Periodic boundary conditions (PBC) are readily applied 
to noninteracting {U — 0) systems whose GS energy is 
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The expression for the infinite chain is shown in Fig. 1; 
E{(1) is the complete elliptic integral of the second kind, 
with q^ = 4/(A2 -f 4). From Eq. (|) the GS electronic 
density on D is: 



71b(A,0)-1 = 



2AX(g) 
7r(4-}-A2)i/2 



(4) 



where K is the complete elliptic integral of the first 
kind. The divergence of dnn/dA at A = signals 
an electronic instability. The behavior of finite rings 
is different and shows 4n, 4n-|-2 effects. The wavevec- 
tor is fc = 0,±27r/A^, ±47r/A^, ...TT. We have energies 
±A at /c = 7r/2 when A^ = 4n and two electrons for 
these orbitals. The degeneracy produces an energy cusp 
at A = 0; no changes discontinuously and the partial 
derivative in Eq. (S) is not defined. Finite rings with 
A^ — 4n-|-2 have nondegenerate GS at A = 0, no cusp 
and finite {dnu/dA)Q. The 4n, 4n-|-2 sequences must 
coincide in the extended chain and do so according to 
Eq. (Q), with continuous no and divergent (9n£)/9A)o. 
Exact U = results illustrate the extrapolation problems 
encountered in interacting chains. 

The full basis of Ho{t, A,U) increases roughly as 4^ 
with A^ and as 3^ when we exclude doubly ionized sites, 
i.e. two electrona-at A sites or two holes at D sites. We 
use VB methodsO to reach A^ = 16 for the full basis and 
A = 22 for the restricted basis. The basis has over 10^ 
singlets peiTr-lO^ Slater determinants with Sz = 0. Exact 
solutior£3S of the extended chain is limited to A = 0, 
the Hubbard model. The GS is a nondegenerate singlet, 
the charge gap is finite for [/ > and there is spin-charge 
separation at large U. Finite t and U always lowers the 
energy in Fig. 1 compared tot = 0. The greatest changes 
occur at A = ±[//2, where t cannot be treated as a small 
parameter. 

In C]\r symmetry, the GS of the interacting systems 
transforms as /c' = tt on the covalent side of 4n rings and 
as fc' = in 4n-|-2 rings. Site energies A > lower the 
symmetry from Cn to Cn/2 and yield a charge-density- 
wave (CDW) GS. The extended system no longer has 
inversion centers between sites, which corresponds to re- 
fiection between sites for finite A^, but retains inversion at 
the sites or, for finite A^, reflection (Ty through the sites. 
With two sites per unit cell, both fc' = and tt transform 
as fc = in the first Brillouin zone. According to reflec- 
tion through sites, the covalent GS of 4n+2 rings is even 
(cr„ = 1, Ai) and that of 4n rings is odd {a^ — —1, A2). 
The GS of 4n rings is degenerate at Ac{U, N), where the 
symmetry switches from A2 to Ai with increasing A, and 
this crossover defines the NIT. In addition to PBC, we 
use antiperiodic boundary conditions (APBC) with re- 
versed sign tiN = — 1 for transfer between 1 and A. This 



corresponds to 4>n+N = — </'n for the on-site wavefunc- 
tions and a periodicity 2N. In terms of VB diagrams, we 
modify a^ for reflection through sites 1 and N/2 + 1 to: 



^.(-1)' 



(5) 



where rii is the occupation number of site 1. The APBC 
GS has ct' = 1 in 4n rings for any A, U. The GS of 
4n-|-2 rings have a crossover from ct' = — 1 at small A to 
cr' = 1 at A > Ac{U,N). The subspaces A[,A2 associ- 
ated with ct' do not coincide with Ai, A2. The paired 
state is unique and even for either PBC or APBC. There 
are two covalent states, the Kekule diagrams for ben- 
zene, with nearest-neighbor pairing of all spins. We de- 
fine \K1) and \K2) as pairing spins at sites 2i— 1, 2i and 
2i, 2i + 1, respectively, for all i. The pairing in \K1) is 
D+A", while the pairing in \K2) is A~ D+. The combi- 
nation \K1) + \K2) transforms as Ai or Aj for PBC and 
APBC, respectively, while the out-of-phase combination 
transforms as A2 or A[. 

The [/ = A = crossover connects electrons paired as 
D+2A~2 oj. DA in Fig. 1. For C/ > 0, the NIT shifts 
to positive U — 2Ac and i 7^ preferentially stabilizes 
the paired GS over the covalent GS because the latter 
has finite probability for adjacent parallel spins that can- 
not transfer under Eq. ([l|). The symmetry changes at 
±Ac{U,N) in rings with either PBC or APBC. Exact 
crossovers are shown in Fig. 2 as C/ — 2Ac(U, N) in the 
C/,A > quadrant for U = 0.5, 1, 2, 3, 4, 5, and 10; 
the inset has U ~ 100, 200, 300 and 00, the last one 
corresponding to the restricted basis. At fixed U and fi- 
nite N, the crossovers are similar for 4n with PBC and 
4n+2 with APBC. The dashed line is an extrapolation to 
the infinite chain discussed below. The covalent region 
is very narrow and the crossovers merge ai U — 0. The 
inset shows that Ac{U, N) is nearly constant for A > 5i. 

We plot Ac{U, N) vs. N~'^ in Fig. 3 and find accurate 
extrapolation at large A, [/. The difference between U = 
300 and the restricted basis is due to small admixtures 
of A^^D^+ at energy U + 2A. The extrapolated limit is 
U — 2Ac = 1-332 in the restricted basis. It has previously 
been estimateda as 1.2-1.3 based on thejSinglet and triplet 
gaps, respectively, of A^ < 10 rings ando at 1.5 based on 
the ionicity up to A^ = 10. Mixing with A^~D^+ grows 
as U decreases, as seen for U = 10, and the functional 
dependence is closer to ~ A^"^ at smaller U ^ 3. Kinetic 
contributions are largest at small U and A, where t is 
comparable to CT energies. The extrapolated (dashed) 
line in Fig. 2 is based on a power law, Ac{U, N) oc N~'^, 
with 1 < 7 < 2 giving the best fit for each U from A^ = 
8 to 16. 

Figures 2 and 3 indicate that, except for A, [/ < 2t, the 
NIT hardly varies with U/A. The relevani_DA systems 
have narrow bands, oxides are modeledQo with wider 
bands, A < t. The restricted basis captures the basic 
physics. We let both A and U diiierge in HQ{t,A,U) 
while keeping F = A — U/2 finiteCI and reference the 
crossover to F — Fc = A — Ac. In the half-filled case, A -^ 



00 ensures an electron at each D site and excludes two 
at any A site. The i = GS has energy -2F, per DA 
on the paired and covalent side, respectively, with F = 
at the NIT. The restricted basis is almost quantitative 
for U > 5, makes larger A^ accessible, and holds at the 
NIT. The related hmit with both C/ and C/ - 2A >> t 
leads instead to a Heisenberg spin chainEJ without charge 
degrees of freedom and does not apply to the NIT. 

The GS expectation value, (n2i_i), for electrons at D 
sites is more accurate than the numerical derivative in 
Eq. (0). Matrix elementsEJ over correlated states can be 
evaluated exactly for finite A^. Results for the restricted 
basis of 4n rings with PBC and 4n+2 rings with APBC 
are shown as a function of F — Fc(A^) in Fig. 4a. The 
crossover generates a jump in no- The charge density 
is continuous in Ai, A2 for 4n rings, or in A'l, A'2 for 
4n+2 rings, but continuing the lines in Fig. 4a through 
the NIT gives an excited-state density. All approaches 
to the NIT described by the hamiltonian in Eq. (|l|) indi- 
cate riD to be continuous when t is finite. As expected, 
the discontinuity in no decreases with N and vanishes 
in the extended chain. The GS of 4n-|-2 rings with PBC 
or 4n rings with APBC remains in Ai or A'-^^, respec- 
tively, for any U, A and the charge density is continuous, 
as shown in Fig. 4b. The NIT defined by the maxi- 
mum of dno/dA is less precise numerically (by ~ 0.02i) 
than a crossover. The curves in Fig. 4b have been ad- 
justed to catch the extrapolation between increasing and 
decreasing series on either side of Ac. Results for the in- 
finite chain are shown as stars that coincide in both pan- 
els. They represent joint extrapolations as either A^^^ 
or N~'^ that give the smaller mean square deviation. In- 
deed, riD is almost quantitatively known from the re- 
quirements that no (22) > 7io(20) on the covalent side 
and nr,{22) < no(20) on the paired side. The present 
estimate for no is 1.314(2) at the NIT, i.e. p = 0.684. 

The restricted basis for U, A >> t fixes n^i = 1.31 at 
the NIT of Eq. (|i|). The slope dno/dA is finite, but 
this is inconclusive by itself. Hiickel rings show similar 
4n, 4n-f 2 behavior and exact A^ = 200 and 400 results are 
indistinguishable from no in Eq. (]^ at the resolution of 
Fig. (4) . The origin must be magnified an order of magni- 
tude to see the divergence of dno/dA — {d'^Ea/dA'^)/N 
at A = 0. This divergence signals the intrinsic instability 
of the U = A — chain to a site-CDW distortion that, 
however, is already broken at finite U = 2Ac in inter- 
acting systems. Extended chains with U > have finite 
dno/dA at the NIT. By contrasl_-the Peierls instability 
to a bond-CDW is unconditionalEHlj for any t/U, because 
dimerization breaks reflection symmetry a^, as experi- 
mentally recognized in the initial TTF-CA studiesB. 

The full basis of the Hamiltonian in Eq. (|l|) is required 
for small U and exact results for no or Ac{U) are lim- 
ited to A^ = 16. We again have discontinuous no{A, U) 
in 4n rings with PBC and 4n-|-2 rings with APBC, and 
continuous no for the opposite boundary conditions. We 
compare in Fig. 5 extrapolated no for U = 2, 5 and 10 
with the exact U = result in Eq. (H). The arrows mark- 



ing the NIT for finite U are based on symmetry crossovers 
and extrapolations similar to Fig. 3. We have increas- 
ing nni^c, U) with U and the limiting value of ^--^1.3 is 
reached hy U = 5. Figure 5 shows the stabilization of co- 
valent states with increasing U and small NIT variations 
for A, C/ > 2. 



III. ENERGY GAPS AND LOCALIZATION AT 
THE NIT 

The excitations of Ho{t, A, U) provide other evidences 
of the NIT. The paired and covalent states are diamag- 
netic and param|affnctic, respectively. A singlet-triplet 
gap, Est, open^ at the NIT between a band insula- 
tor with Est > and a Mott insulator with Est = 0. 
The lowest singlet excitation, M^ss, is between the Ai 
and A2 GS; hence Ess vanishes^ at the NIT. The transi- 
tion is dipole allowed and is formally the CT excitation, 
but Ess rapidly loses oscillator strength on the cova- 
lent side. The charge degeneracy in Fig. 1 at U = 2A, 
i = distinguishes between neutral and ionic complexes. 
The i > gaps near the NIT are not known and their 
simultaneoHS-jOpeaing, as tacitly supposed foc-|apSingle 
transition,!] BtJ EZl is not assured nor agreed onEJ"E3. We 
report exact excitation thresholds near the NIT defined 
by GS crossovers. All symmetry considerations apply to 
the full basis for A > 0. 

Figure 6 reports EssiN) = E2iN) - Ei{N), i.e. the 
energy difference between A2 and Ai GS, in the restricted 
basis as a function of F — rc(iV). Since Ess{N) increases 
in 4n rings for F > Fc and decreases in 4n-|-2 rings, we 
have finite Ess on the paired side. On the covalent side, 
Ess{N) decreases with TV in rings whose GS remains in 
Ai or A[ and increases in rings whose GS is in A2 or 
Aj. Joint extrapolations yield the stars that are consis- 
tent with vanishing Ess in the extended system. Exact 
results to iV = 22 in Fig. 6 are the most stringent limit 
to date, with Ess < 0.05i on the covalent side and finite 
Ess{N) for F - Fc < 0.05i. We note, that at A = 0, 
far on the covalent side, Ovchinnikovt3 found nonpolar 
singlets with zero gap for any U > 0. Far on the paired 
side, we have Ess ~ 2F by inspection. Hence increasing 
A at fixed U in the extended system clearly opens an 
SS gap that is seen to coincide in Fig. 6 with the NIT 
defined by the symmetry crossover. The unconditional 
instability for dimerization on the covalent side is closely 
related to vanishing Ess', the instability is conditional 
for a finite gap. 

The magnetic gap Est is to the lowest triplet for ei- 
ther PBC or APBC. As shown in Fig. 7, Est increases 
rapidly with F > Fc in the restricted basis and is small on 
the covalent side. Open circles represent boundary con- 
ditions with crossovers and systems whose Est increases 
with N at larger F — Fc. Closed circles are for boundary 
conditions without crossovers and show decreasing Est 
with TV. The stars in Fig. 7 are joint extrapolations. 



The bound on Est is Est < O.li for F < Fc and the 
gap opens at Fc or slightly on the covalent side. At fi- 
nite U, the extended system is rigorously known to be 
paramagneticEj at A = 0, with Est = 0, and diamag- 
netic with Est ^ 2A ~ U for A » U . The opening 
of an ST gap with increasing A is assured, and the re- 
sults in Fig. 7 are consistent with Est > at the NIT. 
The concomitant dimerization on the covalent side-fipens 
a magnetic gap, as it is well known in spin chainscl with 
regular or alternating exchanges and triplet spin excitons. 
The charge gap of the Hamiltonian in Eq. d^) is / — A, 
since there is not an explicit Madelung contribution. 
/ — ^ is related to the GS of the cation and anion rad- 
icals, E+{N) and E^{N), respectively, and corresponds 
to charge disproportionation or electron transfer between 
noninteracting systems. 



I{N) - A{N) = E+{N) + E^{N) - 2Eq{N). 



(6) 



At t = and A > 0, we have a paired GS for [7 < 2A 
with I = A-U and A = -A; the fines cross at the NIT 
and the covalent side has / = — A and A = A — U ioi 
U > 2A. For t > 0, the charge gap of free electrons, 
|2A|, follows from the valence and conduction bands in 
Eq. (^); the extended U — system is metallic at A = 
and insulating otherwise. Although not known exactly 
for [/ > 0, the charge gap is readily shown to be large, 
roughly |2A — [/|, far from the-NIT. On the covalent side, 
it becomes the Lieb-Wu gapEd at A = and increases as 
U ioT U > At; on the paired side, all gaps increase as 
2 A — U for A >> U. Finite N leads to charge gaps at 
the NIT in systems wiili discrete energies. 

Both e-h symmetryEj and ab are broken for A > 0, 
but their product reniajns a symmetry operation. Ex- 
tended e-h symmetryEj cuts basis, either full or re- 
stricted, roughly in half and corresponds in the S = 
manifold to the A+ ® B~ and A~ ® B+ subspaces of 
Ho{t,0,U). The GS symmetry does not change at the 
NIT. E-h symm|Gtcy, relates the GS and excited states of 
the radical iona^J'E^. In particular, we have E-(N) = 
E-^-{N) + [/, a general result that holds on adding any 
spin-independent potential to Ho{t, A, U). It follows that 
Eq. (H) reduces to I{N)+A{N) = -f7 for even N and ar- 
bitrary t, A, U. Table 1 reports I{N) - A{N) at U, Ac up 
to A^ = 14 in the full basis and N — 18 in the restricted 
basis. The U — gaps vanish at the crossover, where the 
electron transfer described in Eq. (||) involves degenerate 
orbitals. The charge gaps increase with U but remain 
small at the crossover even for divergent U. The gaps in 
Table 1 follow power laws, N~'>', with 7 < 0.6 and place 
a rough bound of ~ 0.2i on the extended system. The 
charge gap vanishes at most at a single point, Ac{U) or 
Uc{A), that coincides with the NIT within the accuracy 
of finite systems. 

The SS, ST and charge gaps are all finite on the paired 
side, F > Fc. They differ on the covalent side, however, 
where only the charge gap is finite. We associate gaps on 
the paired side with localization. The GS for the |A| >> 



U limit has paired spins on either odd or even sites and 
is manifestly localizedn-Since gapless triplets and singlets 
are firmly establishedEl at A = 0, the GS of iJe(t, A, U) 
for [/, i > is extended at A = and localized at A >> 
U . A localization-delocalization transition between two 
insulators incorporates all aspects of the NIT and the 
vanishing charge gap suggests a metal at the transition. 
We develop these ideas below. 

To show localization on the covalent side, we partition 
Hq in the restricted basis into /iq for transfers between 
sites 2n-l and 2n, as in \K\), and a perturbation V for 
transfers between 2n and 2n+l. We take F = A— [//2 > 
in units of t and solve the 2x2 dimer problem in the 
singlet subspaces of /iq. The exact GS of ho is 



A2 states for PBC. We now have £'o(r, i, /) in units of t. 
The charge stiffness per site is Xcs = {d'^Eo/dp)o/N . It 
is finite in conductors and vanishes in insulators. At the 
crossover, the proper zcroth-order GS of 4n rings is the 
odd linear combination of the Ai and A2 GS and 



Xcs (re) = 



.(1) 



N/2 



|Go(r)) = n [4-i.««: 



2i-l,/3 



COS( 



V2 



'2'i-l,a"2i,/3 



^2f 



(7) 



where |0) is the vacuum state and tan 20 ~ \/2/r governs 
the mixing of \DA) and the singlet linear combination of 
\D+A-). The zeroth-order energy is -T -~ (F^ + 2)^/2 
per dimer. The opposite choice of 2n, 2n+l for dimers 
has the same energy but admixes \A~D^) singlets, as in 
\K2). Each dimer has a triplet with excitation energy 
ET = r + (r^ + 2)1/^ a singlet at es = 2{r'^ + 2f''^ and 
strictly confined electrons. The perturbation 



V = - 






{o.2i^a"'2i+l,cr + <Ai+l,a'^'2i,<y) 



(8) 



is necessarily small when T is large. To second-order in 
y, the energy per dimer is 



.(0) 



e(2) = -r-(r2 + 2)1/2- 



cos 



(r2+ 2)1/2 



(9) 



Electrons are now confined to adjacent dimers that are 
connected by virtual excitations. As shown in Table 2, 
Eq. (P) is nearly quantitative as close to the NIT as F = 2. 
Localization to adjacent dimers approximates the exact 
solution of Eq. (|l|), which for A^ = 22 is a linear com- 
binations of over 10^ singlets. Rapid convergence with 
N also points to localization on the paired side and is 
seen for tid in Fig. 4, -Ess in Fig- 6 and Est in Fig. 
7. Successive orders in V increase by one the number of 
coupled dimers. Such an expansion fails at the NIT or 
on the covalent side. 

To see if the system is metallic at the transition, we 
compute the charge stiffnes£3 relevant to the Hamilto- 
nian in Eq. (m)-^!!!^ property has been applied to inter- 
acting fermiona23c3 in one dimension. The perturbation 
is a phase factor exp (±i/) in Eq. (IT]) for transfers to the 
right and leftcj. 



V{f) = {l-cosf)v++iv^sinf 



(10) 



The first term of Eq. (IT]) is —v+^ while v- has oppositely 
signed transfers to the right and left and connects Ai and 



Ndp 



,(2) 



= p« 



.p(2) = 



1^ 

N 



F,(n 



.(2) 



(11) 



Here -n}^' and n\j' are the electron densities at donor sites 
in the Ai and A2 subspaces, respectively, and p^^> and 
p(2) are the corresponding bond orders, with 2p defined 
as the GS expectation value of v+ , which we evaluate in 
the restricted ba,sis. The donor densities at the NIT are 
-^I^fst^juirbF^P 4 a|fl^ have opposite N dependence in Ai 
and A2. The value of F^ = U/2 - A^ = -0.666 follows 
from Fig. 3. We have a poor metal: XcsiXc) ^ 0.74 
is ~60% of A/tt value for free electrons aX A — U — 
in Eq. (Q). For F 7^ IL, second-order perturbation 
theory in / becomes exactE3. Fig. 8 shows that Xci(r) is 
exponeiitially small on either side of the NIT. The charge 
stiffnesa^a of the Hubbard model has a similar peak at 
[/ = that narrows with A^ and becomes a (5-function 
in the infinite chain. Hence we expect Xcs(r) to vanish 
except at Tc in the extended system. 

A metallic point connecting inS|»iating phases has been 
recently discussed for i?o(ij A, [7)0 and for the following 
half-filled system of spinless fermionso, 

H = - ^(a^+ai+i + a^jfli) + '^{VniUi+i + Wnini+2) 

i i 

(12) 

Large V > favors a GS without adjacent occupied sites, 
while large W favors one with adjacent filled and empty 
sites along the chain. The t — Q crossover occurs at V = 
2yV, with adjacent electrons and holes for V < 2W that 
resemble D and A sites, respectively. The V > 2W GS 
has alternating filled and empty sites that, taken in pairs, 
correspond to D"*" and A~. The crossover is not precisely 
at V = 2yV, presumably due to different bond orders in 

the two GS. It shifts to W - V/2 0.6 in units of i, 

very close to T^- Transfers between two sites in Eq. (|l^ ) 
differ from spin degeneracy in Eq. (nl), however, aod the 
models do not map into each other. Exact resultsa to A^ 
— 40 for Eq. (Jl^) are comparable to A^ = 20 for the full 
basis of Eq. (|l^ The striking similarity between Fig. 8 
and the charge stiffness of Eq. (na) up to A^ = 40 suggests 
a similar interpretation. 

IV. DISCUSSION 

The modified Hubbard model in Eq. (|l|) has many ap- 
plications to both theory and experiment. With vari- 
ations, it is suitable for modeling valence transitions, 



2) 



excitation thresholds, electronic or structural instabili- 
ties, among other topics. Its two parameters, U /t and 
A/i, encompass the Hubbard model (A = 0) at half or 
other filling, two bands at C/ = and localized dimers for 
A >> C/. At fixed U and t, increasing A > generates 
a neutral-ionic transition whose characterization is the 
principal goal of this paper. The NIT of the Hamiltonian 
in Eq. (|l|) has continuous ionicity given by Eq. (g) and 
excitation gaps for singlets, triplets and charges that are 
not known exactly. Previous approximations have been 
developed separately for no, Ess, Est, the charge gap, 
the GS at the NIT, instabilities, etc. Our collective anal- 
ysis of symmetry crossovers, excitation thresholds and 
GS properties incorporates computational advances and 
yields better estimates for extended systems. 

Finite-size results require extrapolations whose ac- 
curacy improves with N. We followed the NIT of 
i/o(i, A, U) through the symmetry crossover of the GS, 
the charge density ud, the excitations Ess, Est, and 
the charge gap and stiffness. Larger N is accessible in 
the restricted basis allowing for an accurate estimate of 
Tc == — 0.666i from the crossover and setting stringent 
limits of ^ O.lt for the opening of all three gaps at this 
position. Thus the numerical results point to a single 
transition. The NIT of the modified Hubbard model is 
continuous, as previously found, and is marked by the 
opening of singlet, triplet and charge gap on the paired 
side. There is no gap in the singlet or triplet manifold 
on the covalent side. The interacting system is known to 
have a delocalized GS at A = 0, the Hubbard limit, and 
localized GS for A >> U, the paired limit. We identify 
the NIT at ±Ac([/, t) as the appearance of a localized 
GS. ^ 

Resta and Sorella discusa^ polarization and metallic 
behavior at the NIT in the context of oxides, with tg = 3.5 
eV, A' — 2.0 eV and variable U in Eq. (|l|). The crossover 
in A^ = 8 rings, at U/to = 2.27, is used to estimate the 
polarization of extended systems. Since A' /to = 0.571 
corresponds to large t, the crossover is near the origin 
of the A, U plane in Fig. 2 and there are substantial 
finite-size effects. The iV = 8 result in Fig. 2-4'ields 
Uc — 2.27, in quantitative agreement with ref.cfl, but 
larger N up to 16 extrapolate to larger Uc/t = 2.70 for the 
extended system. Such corrections are consistent with 
A^ ~ 10 results on CT complexes. The GS polarizability 
is a new approach, different from the charge stiffness, to 
the identification of metallic behavior. 

The GS density is no = 1.314 at the NIT of the re- 
stricted basis, when one electron is always confined to D. 
The spin degeneracy of D+ or A^ spoils exact analysis. 
The degeneracy of charge and spin excitations at the NIT 
gives a simple, heuristic interpretation: nu = 4/3 is the 
result for equal weights of molecules and spin-1/2 radical 
ions. Equal weights at the NIT can be justified rigorously 
at [/ = for electrons or for spinless fermions, but not 
in the restricted basis. We found Tc = -0.666 in the re- 
stricted basis and use this value in |Go(r)), the dimer GS 
in Eq. (0); the paired-state amplitude is cos^ </)= 0.287, 



which corresponds to n^i — 1.287. Dimers capture most 
of the configuration mixing of the extended system. The 
full basis has contributions from D^+ and A^^ diagrams, 
which as seen in Fig. 5 reduce no compared to the re- 
stricted basis. 

Peierls-Hubbard models are widely applied to struc- 
tural instabilities. The stability of the GS to a pertur- 
bation can be formulated in terms of susceptibilities, x, 
that are formally given by the exact eigenstates \F) and 
energies Ep of the hamiltonian in Eq. (|l|) . The perturba- 
tion is written as the product 9Q, where 9 is the relevant 
operator for coupling to Q, and the corresponding x is 



X« 






2E 



\{G\0\F)\' 
Ef-Eg 



(13) 



Since the sum is over the excited states of the unper- 
turbed system, the eigenstates of the uniform chain in 
Eq. (H) suffice for the stability of the modified Hubbard 
model. The charge stiffness in Eq. dl3) is x ^^ ''*5- 
spect to a magnetic field perpendicular to the ringLJ and 
gives information about current fiow. Structural tran- 
sitions are investigated by introducing phonons as Q- 
perturbation. The Peierls instability for dimerization in- 
volves fc = phonons with 9 representing the staggered 
bond-order operator (the first term in Eq. (U), augmented 
by a (—1)* factor). This operator breaks inversion sym- 
metry at the sites and mixes Ai and A2 singletaj. Van- 
ishing Ess on the covalent side of the NIT then indi- 
cates a divergent x and the unconditional instability of 
a lattice with harmonic potentialaa. On-site (Holstein) 
phonons couple instead to CDW operator no- Since 
XA — drto/dA is finite at the NIT, except for A = 0, the 
corresponding instability is conditionals; the NIT marks 
the maximum XA, i-e. the maximum dnu/dA, as dis- 
cussed under Fig. 4. 

We turn next to open or controversial aspects ot-tha 
NIT of the modified Hubbard model. Some authorg2a^E2l 
have proposed two transitions related to the closing 
of charge and spin gaps, respectively; a spontaneously 
dimerized phase then separates a band insulator corre- 
sponding to Ihe paired GS and Mptt insulator on the 
covalent sidecj. The suggestion£j for another transi- 
tion rest on the analogy with spin-1/2 Heisenberg anti- 
ferromagnetic chains with frustration due to a second- 
neighbor exchange J2- The Kekule diagram \K1) or 
|Ar2) is the-|exact GS at J2 = ^i/2, as recognized by 
Majumdarcil. There is no exact mapping of Eq. (1) 
into such a spin chain, not even at large U, but the 
GS of related models with, for example, second-neighbor 
transfers have not been studied in detail. Our exact re- 
sults for iJo(t, A,[/) with finite A show that the max- 
imum of dp/dA, the closing of the singlet and triplet 
gaps, and the vanishing of the charge gap coincide at 
the NIT within ^ O.lt. Finite systems cannot specify 
transitions, but provide some constraints. The Mott in- 
sulator dimerizes spontaneously,for any t/C/, as discussed 
for the spin-Peierls instabilityc3 of Heisenberg antiferro- 



magnetic chains. The dimerization ampUtude becomes 
very small for U >> t and J = t'^/U, since the electronic 
stabilization is less than J, but the singularity actually 
increasesCJ; the GS energy in Eq. (||) at t/ = A = goes 
as S^lnS for alternating t{l ± S) along the chain, while 
the GS of the spin chain with alternating J(l ± 6) goes 
as iJ^/"^ In (5cJ. Such considerations apply to Eq. (|^) in 
the covalfijit limit T = U — 2 A >> t where, as noted 
originallyaa, we have a Heisenberg chain with J — t^ /T. 

The charge gap is a recejat— topic and is expected to 
have a minimum at NITEj'LJ-c3. We find finite minima 
in interacting systems with finite N. As already noted, 
the polarizabilitycJ and the charge stiffness in Eq. ( p^ ) 
and Fig. 8 give independent indications of a metallic 
GS at the NIT. We consequently expect vanishing charge 
gap there. The metal separating two insulating phases 
is extremely fragile: not only is it restricted to a single 
point Ac(C/), but it is unconditionally unstable to dimer- 
ization. Moreover, extending the model in Eq. (yj) to 
include intersite e-e interactions or on-site phonons pro- 
duces a discontinuous NIT above some critical coupling, 
which excludes a metallic phase even for a rigid lattice. 
By contrast, a metallic GS at the NIT of Eq. (^ is fairly 
robust. We have a simple half-filled band at f7 = and 
a correlated metal persists to arbitrarily large U. At 
the NIT, A counterbalances U: the charge distributions 
DA and D+A^ are almost degenerate and hence strongly 
mixed by any finite t. A metallic state at NIT is not con- 
sistent with finite triplet excitation there. 

The degeneracy of charge and spin excitations is char- 
acteristic of the NIT and appears already in the t — 
limit of Fig. 1. Finite A completely spoils spin-charge 
separation at the NIT of Eq. (|l|). Vanishing Est is closely 
linked to magnetic susceptibility of Hubbard or Heisen- 
berg chains. Since a singlet can always be constructed 
from two triplets, vanishing Ess follows immediately and 
is associated with even-parity spin waves in A = sys- 
tems with e-h symmetry. This symmetry is broken in 
Eq. (P or its extensions, and the CT excitation connect- 
ing the Ai and A2 GS is dipole allowed. Spin-charge sep- 
aration is regained on the covalent side when the charge 
gap exceeds a few t, much as in Hubbard models for 
U > t: Exact separation requires infinite U, but U > At 
suffices in practice. 

To summarize, we have extended exact solutions of 
the modified Hubbard model in Eq. (|]) to larger sys- 
tems, identified the NIT with symmetry crossovers in 
rings with periodic or antiperiodic boundary conditions, 
and found the charge density, excitation thresholds and 
susceptibilities at the NIT. We find a continuous NIT and 
tighten considerably the extrapolated limits for the infi- 
nite chain. Our results indicate a T=OK transition with 
vanishing singlet and triplet gaps on the covalent side, 
vanishing charge gap and metallic GS at the NIT, and 
finite singlet, triplet and charge gaps on the paired side, 
whose localized GS is confirmed. We associate the NIT of 
the model in Eq. (|l|) with a transition from a delocalized 
(small A) to a localized (large A) GS. Accurate analysis 



of the hamiltonian in Eq. (|l|) is required to model valence 
transition in charge-transfer salts or metal oxides where 
long-range Coulomb interactions and e-ph coupling have 
to be considered explicitly. 
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TABLE L Exact charge gap, I ~ A in Eq. (|6|), at the neu- 
tral-ionic transition of the model in Eq. (1), for rings of A'' 
sites, with t = 1, variable U and Ac{N) at the crossover of 
4n (4n+2) rings with periodic (antiperiodic) boundary condi- 
tions. 



N 


[7 = 2 


10 


cx) (restricted basis) 


8 


0.1893 


0.7162 


0.8486 


10 


0.1803 


0.6548 


0.7614 


12 


0.1718 


0.6010 


0.6909 


14 


0.1614 


0.5555 


0.6342 


16 






0.5874 


18 






0.5469 



TABLE H. Approximate ground states energy per dimer, 
Eq. (9), of the infinite chain and exact results for Eq. (1) 
with iV = 12, i = 1, r = A - Ac(12, U) and U = 10 and oo 
(the restricted basis). 

F U = 10, exact U = oo, exact i7 = oo , Eq. (9) 
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FIG. 3. Size dependence of the gs crossover between A'^ — 
8 and 14 at U/t = 3, 10, 300 for the fuU basis of Eq. (1) and 
up to A^ = 20 in the restricted basis with infinite U. 



FIG. 1. Ground-state energy per site, E/N, of the modi- 
fied Hubbard model, Eq. (1), as a function of the site energy 
A for free {U = 0) or interacting {U > 0) electrons in the limit 
of no overlap (f = 0), with valence transitions at A = ±C//2 
in donor-acceptor stacks. The t = 1 curve for free electrons is 
Eq. (3). 




FIG. 2. Ground state crossovers, [/(Ac, A'^), of A'^-site mod- 
ified Hubbard rings (Eq. (1)) with periodic and antiperiodic 
boundary conditions, respectively, for A' = 4n and 4n-f 2. The 
dashed lines are N ^ oo extrapolations discussed in the text. 
The inset shows the large-A behavior and the restricted basis 
at A ^ oo. 
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FIG. 4. Ground-state electron density, ud, of A'-sitc mod- 
ified Hubbard rings, Eq. (1), in the restricted basis. The 
boundary conditions in (a) produce symmetry crossover at 
rc(A'^), the vertical line, where no increases discontinuously 
with r and the smallest jump occurs for A' = 20. The bound- 
ary conditions in (b) with the same rc(A'') result in continuous 
no with increasing dno/dV up to A' = 22. The stars are joint 
N ~* oo extrapolations of (a) and (b) discussed in the text. 
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FIG. 5. Ground-state electron density, no, of modified 
Hubbard models, Eq. (I), with V jt = 0, 2, 5 and 10. The 
exact [7 = result is Eq. (4); [/ > points are N -^ oo 
extrapolation of no based on the full basis up to A'' = 16; the 
arrows mark the neutral-ionic transition found as in Fig. 3. 



FIG. 7. The singlet-triplet gap. Est, near the NIT of Eq. 
(1) up to Af = 18 in the restricted basis. Open and closed cir- 
cles refer to boundary conditions with and without crossovers, 
respectively, and stars are joint N ^ oo extrapolations based 
on both. 
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FIG. 6. The singlet-singlet gap. Ess, near the NIT of Eq. 
(1) up to A'^ = 22 in the restricted basis. Boundary con- 
ditions leading to crossovers are shown as open circles and 
\Ess\ is the excitation for F < Fc. Boundary conditions with- 
out crossovers are shown as closed circles. The stars are joint 
N —> oo extrapolations based on both. 



FIG. 8. Charge stiffness, Xcs(r) in Eq. (12) and as dis- 
cussed in the text, near the NIT of N-site model, Eq. (1), 
in the restricted basis; Xcs(Fc) is ~ 60 % of the free-fermion 
value. 
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